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Abstract 

Suppose we have a Nikishin system of p measures with the kth generating measure of 
the Nikishin system supported on an interval Ak C R with Ak PI A^+i — for aU k. It 
is weU known that the corresponding staircase sequence of multiple orthogonal polynomials 
satisfies a (p + 2)-term recurrence relation whose recurrence coefficients, under appropriate 
assumptions on the generating measures, have periodic limits of period p. (The limit values 
depend only on the positions of the intervals A^.) Taking these periodic limit values as 
the coefficients of a new (p + 2)-term recurrence relation, we construct a canonical sequence 
of monic polynomials {P„}^oi the so-called Chebyshev-Nikishin polynomials. We show 
that the polynomials P„ themselves form a sequence of multiple orthogonal polynomials 
with respect to some Nikishin system of measures, with the feth generating measure being 
absolutely continuous on Afe. In this way we generalize a result of the third author and 
Rocha [52] for the case p = 2. The proof uses the connection with block Toeplitz matrices, 
and with a certain Riemann surface of genus zero. We also obtain strong asymptotics and an 
exact Widom-type formula for the second kind functions of the Nikishin system for {P„}5^q. 

Keywords: Multiple orthogonal polynomial, Nikishin system, block Toeplitz matrix, 
Hermite-Pade approximant, strong asymptotics, ratio asymptotics. 
2010 AMS classification: Primary 42C05; Secondary 41A21. 

1 Introduction and statement of results 
1.1 Nikishin system 

Let p e Z>o and let Ai, . . . , Ap C M be compact intervals such that 

AfcnAfe+1^0, fc = l,...,p-l. 

Assume that for each k £ {1, . . . crfc is a finite positive measure supported on A^ w\t\v density 
(tJ.(x) > for a.e. x G Afc. We denote w\t\v 

= (pii,...,pip) =A/'(cri,...,crp) (1.1) 
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the Nikishin system generated by the measures ak- The construction of such a system is based 
on the following "product operation" for measures; given measures CTq, ct/j, supported on disjoint 
compact intervals on M, set 



d(crQ,cr/3)(a;) 



X — t 



This defines a new measure whose support coincides with that of CTq. The system M in (11.11) is 
then defined as follows (the notation was introduced in [15 ): 



^i:=o-i, A*2 := (o-i,cr2), A^s := (o-i, 0-2, ^3) = (di, (ct2, CT3)), 



:= (cti, (0-2, • • • , CTp)). 



Thus the measures ni, . . . , Hp are all of fixed sign and supported on Ai. Nikishin systems were 
introduced in [24] . These systems (and variations of them) have attracted an ever increasing 
interest during the last decades, due to their many theoretical and practical applications, see 

e.g., [iiiiiaiaiaiHiiniiiiiiiiiiiiisiiiniiiniiiiiEaiMiES]. 

For n G Z>o we define the multi-index 



(rii, . . . , rip) := (to + 1, . . . , m + 1, TO, . . . , m) e Z' 

k times p — k times 



P 

>0' 



(1.2) 



where to G Z>o and k G {0,1, ... ,p — 1} are such that n = mp + k. Note that we have 
\n\ :— m + n2 + ■ ■ ■ + Up = n. 

Let (Qn)5?Lo be the diagonal sequence of multiple orthogonal polynomials associated with 
the Nikishin system Ai. That is, the polynomial Q„ is the monic polynomial of degree n that 
satisfies for each k = 1, . . . ,p the orthogonality conditions 



Ai 



i(a;) d^fe(x) = 0, / = 0, 1, . . . , rife - 1, 



(1.3) 



where is the fc-th component of the multi-index n in (jl.2p . The existence and uniqueness 
of the sequence (Qn)^o follows from the (weak) perfectness of Nikishin systems, see e.g. [25]. 
(See also [13] where it is shown that Nikishin systems are perfect.) It is well known that the 
polynomials Q„ satisfy a {p + 2)-term recurrence relation 

(z), (1.4) 

with initial conditions 

Q-p = Q~p+i = ■ ■ ■ = Q-i = 0, Qo = Ij 

and with an,n-p 7^ for all n > p. The recurrence coefficients a„^m in (|1.4[) can be viewed as the 
entries of the banded Hessenberg operator 



A 



ai,o 0,1,1 



flp+l,! 1p+l,2 



V 



(1.5) 



Thus A is a semi-infinite matrix with one superdiagonal, filled with I's, and with p subdiagonals. 
All the other entries of A are equal to zero. 
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1.2 Chebyshev-Nikishin polynomials 

It was proved in [5 that the fohowing ratio asymptotic formulas hold: 

Qmp-^k (-2' ) 



l...,p, 



(1.6) 



where the limits are valid miiformly on compact subsets of C \ Ai. The limiting ratios Fk{z) will 
be described in more detail in Section H) 

The ratio asymptotics (jl.6p imply, see [2|, that for each fixed {i, j), < j < p—1, j < i < j+P 
the following limit exists: 



lim a. 



(1.7) 



For all m G Z>o, we take b. 



mp+i.mp-\-j 



bij. This means that the banded Hessenberg operator 



A in (jl.5l) is a compact perturbation of the tridiagonal block Toeplitz operator 



T = 



Bi Bq B^i 

Bi Bq B^i 



El Bq 



\ 



(1.8) 



where the blocks B^., are of size p x p and given by 



Bo 



1 



\ 



\bp- 



/bp,( 



1,0 "p-1,1 



bp-i^p-ij 



■J 



bp,i 
bp+1,1 



bp.p—i 
bp+i^p-i 

b2p-i,p-i I 



(1.9) 



/O ... 0\ 



B-i 





0/ 





VI 

We view the diagonals of T as infinite periodic sequences with period p. 
Theorem 1.1. The following relations hold for all k,j and all I G {1, . . . ,p — 1}: 

bk,k) 



(1.10) 



bk+i,k-i — bk,k-i-i = bk,k-i{bk-Lk-i 



(1.11) 

(1.12) 



Theorem 11.11 generalizes a result in (55] for the case p = 2. It will be proved in Section [T] 
However, Theorem ll.il will not be used in the rest of the paper. 

We define the Chebyshev-Nikishin polynomials [22^ as the monic polynomials {Pn)^=o ^-sso- 
ciated with the operator T in (|1.8p : that is. 



F„(z) =det(z/„-r„). 



(1.13) 



where /„ denotes the identity matrix of size n, and T„ is the n x n principal submatrix of T. 
The polynomials Pn satisfy the recurrence relation 



zPn{z) = P„+l{z) + 6„,„P„(z) H h bn^n-pPn-p{z) 
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with the initial conditions 



P-p = P-p+i = ■ • • = -P-i =0, Pq = 1- 

Note that for p = 1 and Ai = [—2,2], the polynomials P„ reduce to the classical Chebyshev 
polynomials of the second kind for this interval. This explains the name Chebyshev-Nikishin 
polynomials |22) . 

1.3 Nikishin system generated by the Chebyshev-Nikishin polynomials 

It is easy to see that there exist positive measures {i>i, . . . , Vp) supported on Ai so that 

-1^=/ -^di^kit), k = l,...,p. (1.14) 

Fk{z) Jai z-t 

This is a consequence of the ratio asymptotics in (II. 6p and the interlacing of the zeros of Qmp+k (z) 
and Qmp+k~i{z) Alternatively, the integral representation (jl.l4l) can be obtained from the 
analytic properties of "p^^ using [18l Theorem A. 6] (see (j4.2l) below) which ensures that the 
measures i^k^k = 1, . . . ,p, are absolutely continuous with respect to the Lebesgue measure. 

The main goal of this paper is to establish the following two theorems on the Chebyshev- 
Nikishin polynomials Pn{z). 

Theorem 1.2. Let {vi, . . . ,Vp) be the measures on Ai defined in (jl.l4p . Then the Chebyshev- 
Nikishin polynomials Pn satisfy the orthogonality conditions 

/ Pn{x)x'- dvk{x) =Q, Z = 0, . . . - 1, fc = l,...,p, (1.15) 
where Uk is the k-th component of the multi-index n in (|1.2D . 

Theorem 1.3. The measures {i>i, . . . ^Vp) form a Nikishin system on (Ai,...,Ap). That is, 
there exist measures pk supported on A^, k — 1, . . . ,p, such that (z^i, . . . ^Vp) — J\f{pi, . . . , Pp). 
The measure pk is absolutely continuous on Ak, with its density p'^ix) being non-vanishing and 
of a fixed sign ( either positive or negative ) on the interior of A^ ■ 

Theorems 11.21 and 11.31 will be proved in Sections [5] and [S] respectively. For p = 2 these 
theorems reduce to a result in [52]. Our proofs will be markedly different from the ones in |22| . 
A generalization of Theorem 11.21 to other systems of multi- indices n will be given in Section [51 

1.4 About the proof of Theorem 11.2} ratio asymptotics 

Let 

F{z,x) := -B^i+ Bo + zBi~ xip (1.16) 

be the block Toeplitz symbol associated with T in (|1.8p . and let 

f{z,x) ■.^detF{z,x). (1.17) 

The algebraic equation f{z, a;) = has p + I roots zi{x), . . . , 2:^4-1(2:) (counting multiplicities) 
that we label so that 

\ziix)\ < \z2{x)\ < ■ ■ ■ < \zp+i{x)\, xeC. (1.18) 
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We introduce the sets 

Tk:={x£C:\zk{x)\ = \zk+iix)\}, k = l,...,p. (1.19) 
In the proof of Theorem 1 1.21 the following result will play a pivotal role. 
Proposition 1.4. We have that 

Tfc^Afc, k = l,...,p, (1.20) 

and 

j.^^ Prnp+kix) ^ ^.^ .Rjh^+jM. ._F,{x), fc-l,...,P, (1.21) 

m^oo Frnp+k-l[.x) ™-^oo C^mp+fc-l (x) 

uniformly on each compact subset of C\ (Ai U A), with A a set of finite cardinality. Here we 
use the notation in (|1.6p . 

Proposition 11.41 will be proved in Section The equality (|1.20p was already obtained by 
Aptekarev [1] Prop. 2.1-2.2] but we will provide an alternative proof. 

Remark 1.5. Once we have Theorems 11.21 and 11.31 to our disposal, it will follow that the poly- 
nomials Pn are of multiple orthogonality with respect to a Nikishin system. Equation (|1.2ip is 
then a consequence of the main result in [5, . It also follows that the zeros of P„ are simple and 
lie in the interior of Ai, with the zeros of two consecutive polynomials interlacing [S] Theorem 
2.1]. Therefore in (|1.2ip the convergence is uniform on each compact subset of C \ Ai. 



1.5 Strong asymptotics, Widom-type formulas, and the second kind 
functions 

As a consequence of Widom's determinant identity for block Toeplitz matrices [27l Section 6] , 
we have the following exact formula for Pmp+k{x) (see [11, Section 4]) 

for all m e Z>o and fc = 0, 1, . . . ,p — 1, where we write 

p-i 

/p = n (1-23) 
i=0 

and where we denote with fl^'-'^+^l the submatrix of F in (|1.16p obtained by deleting row p and 
column A; + 1. Note that we label the rows and columns of F as 1, 2, . . . i.e., we start counting 
from 1 rather than 0. In writing (jl.22p . we are implicitly using the fact that fp ^ 0. This will 
be justified in Section |31 see (|4.5|) . 

Equation (jl.22p holds for all x for which the roots Zj{x) are pairwise distinct. If a; £ C is 
such that two roots Zj(x) and Zk{x) with j ^ k are equal (there are only finitely many such x) 
then (jl.22l) remains true provided we replace the right hand side by its limiting value. 

As a consequence of (|1.22p we have the strong asymptotic formula 
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for any fixed fc e {0, . . . ,p — 1}, uniformly on compact subsets of C \ Fi Section 4]. 

In this paper we will obtain similar Widom-type formulas for the second kind functions of the 
Nikishin system generated by the P„. Recall that the second kind functions '^njiz) are defined 
recursively by |15j 

^n,oiz) ■■= Pn{z) (1-25) 

and 

*n,K^)= / "^"-''f^ d^t), (1.26) 

for I = where pi is the measure supported on A; described in Theorem 11.31 The 

statement of the Widom-type formulas for requires some extra notations and is deferred 

to Section [51 



1.6 Outline of the paper 

This paper is organized as follows. In Section [2] we recall some general results on multiple or- 
thogonality relations for banded Hessenberg operators. In Section [3] we apply these results to 
the Che byshev- Nikishin polynomials P„. Section 2] proves Proposition 11.41 on the ratio asymp- 
totics of the polynomials P„. In Section [5] we use these considerations to prove Theorem 1 1.2 1 In 
Section E] we prove Theorem 11.31 In Section [7] we prove Theorem 11.11 on the relations between 
the coefficients bij in (II. 7p . In Section |8] we obtain strong asymptotics and exact Widom-type 
formulas for the second kind functions \['„^;(z) in (|1.25|) - (|1.26p . Finally, Section [S] discusses the 
generalization of our results to some alternative systems of multi-indices n. 



2 Multiple orthogonality relations for banded Hessenberg 
operators 

In this section we recall some general results on multiple orthogonality relations for banded 
Hessenberg matrices, following Kaliaguine [1^], see also Van Iseghem [22 and Kaliaguine [T7]. We 
consider a banded Hessenberg matrix A withp-|-2 diagonals of the form (|1.5p . We associate to A 
the following sequences of polynomials that we denote by (g„)„ and (pn^)n, where j G {1, . . . ,p}. 
We require these polynomials to satisfy the recurrence relation 



with initial conditions 



n > 0, 



n = 




-p+l 


-p + 2 .. 


. -1 





Pn — 


1 








. 





(2) 
Pn = 





1 





. 





Pn — 








1 


. 





Jp) 1 

Pn — 











. 1 





Qn = 











. 


1 



(2.1) 



(2.2) 



In our definition, we take Oi.j — —1 for < z < p — 1 and j < —1. Observe that for n > 0, is 
a polynomial of degree n and pn^ is a polynomial of degree n ~ 1 for all j G {1, . . . 
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Let An be the truncation of A to the first n rows and columns. We define 

D„(z) :=det(z/„ -AO, (2-3) 
Di^\z) det(z/„_, - {Anf-''''-'^), 1 < j < P, (2.4) 

where Ik denotes the identity matrix of size k, and Af t^' - '-''^'- -'^l denotes the submatrix of M 
obtained after deleting the first j rows and columns. Of course, we set 

Do{z) ^ d['\z) EE D^^\z) ^ ... EE D^;\z) ^ 1, 

and we take Dn{z) = 0, n < 0, and Dn ^z) = 0, n < j. 
Expanding D„ by its last row, it is easy to check that 

Qn = -D„, n > 

since both sequences satisfy the same recurrence relation and initial conditions. We also have 
Lemma 2.1. The following relation holds for every j G {I, . . . ,p}: 

piP = Di'^ + + . . , + jjU) ^ > 0. (2.5) 

Proof. It easily follows from the definition (|2.4I) that the sequence of polynomials Dn"^ satisfies 
the recurrence relation 

zD^J^z) = dI^I,{z) + a„,„ Dl^\z) + a„,„_i Dl^l^iz) + ■■■ + a^^^-p D^^l^iz), n > j. (2.6) 

Since p[^\z) = 1 = d[^\z), and pi^^ = oiP = for 1 - p < n < 0, it follows from (EH]) and 
SIM that p'n^ = dIP for aU n > 0. 
Next we show that 

Pii'' (z) - pI^"^^ {z) = D^i^ [z) for aU n > 0. (2.7) 

Evidently, 

zij^ ~P^r'^) = T^ii-Pi;^ +<^n,n{j^ -P'i''^) + - ■ •+a„,„_p(pi^2,-p!^r;^ n > 0. (2.8) 

Recall that the initial conditions (|2.2p hold, in particular P^f'li-p = 1 and P^-'^2 -p = 1- -l^y 
definition, afe.j = —1 for < A: < p — 1 and j < —1, so if we set n = in (12. 8p we get 
Pi —Pi — 0. If we continue evaluating n from 1 to j — 2 in (j2.8p . we obtain p„ — = 

for 2 < n < j - 1. Therefore (EJ]) holds for < n < j - 1. 

If we now put ?i = j - 1 in we get p'j'^ ~ p'-^^^^ = 1 = o'fK And now (ET|) will follow 

immediately from (|2.8p , (j2.6p , and the fact that pn ~ Pn 

= for j — p < n < j . Therefore 
(|2.5p can be proved now by induction on j. □ 

We recall some results due to Kaliaguine jT^l HZ] , see also Van Iseghem [5^ . In order to state 
these results, we introduce some notation. We will assume that the coefficients of the matrix A 
are uniformly bounded. Let {e„}^Q denote the standard basis in P, and consider the following 
resolvent functions 

5,(z) :=(i?,ej_i,eo), 1 < J < P, (2.9) 
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where Rz = {zl — A) ^ is the resolvent operator, and (•, •) is the standard inner product in P . 
We also define 

0,(z) :=5i(z)+.g2(z) + ■■•+. gj(^), 1 < J < p. (2.10) 
For each j = 1, . . . ,p, we introduce a linear functional Lj defined on the space of polynomials by 

L,(z") = (A"v„eo), 



of "moments" for (pj, in the sense that 



where Vj := eo H h ^j-i- Observe that the sequence (Mi,n)5^Lo = (^j(^"))5^Lo is the sequence 

oo 



n=0 

for all z G C sufficiently large. 

Theorem 2.2 (See [Tni[2S])- a) The vector of rational functions 

(1) (2) (p). 



/ Prt Pn Ph \ 



is a Hermite-Pade approximant to the system {(f>i, 02, . . . , (/>p), with respect to the multi-index n 
in p.2p . r/iai is, /or eac/i j = 1, . . . ,p, we have 



Qn{z)<l>j{z)-p^^\z)^0(^-^y z 



— ?> oo. 



where nj is the j-th component of n. 

b) For each j G {1, . . . the polynomial q„ satisfies the following orthogonality conditions: 

Lj((7„z')=0, l = 0,...,nj - I. 

Remark 2.3. Making use of Theorem 12.21 and Lemma I^TTl by linearity we obtain that 



^ T^n R^n R^n 

is a Hermite-Pade approximant to the system of functions 

(51,51 + Ci,i52,5l + Ci,252 + 22,233, ■■■,91+ Ci,p_i52 H h Cp-i^p-igp), 

with respect to the multi-index (jl.2l) . Here, c,;j denote arbitrary constants. 

The recent survey 4 contains information on the latest developments in the theory of 
Hermite-Pade approximation, multiple orthogonal polynomials, and their applications to ran- 
dom matrix theory. 



3 Chebyshev-Nikishin polynomials 

In this section we apply the results of Section [2] to the Chebyshev-Nikishin polynomials {Rn)'^^Q 
defined by (ll.7p - (|1.13p . To this end it is convenient to use an alternative determinantal formula 
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for P„. In fact, if we applied the results of Section [2] directly to the operator T in (|1.8I) then we 
would need to control the limiting ratios 



det (z/„ 



{Tn) 



for n — oo and for fixed j ~ 1, . . . ,p. These limiting ratios are hard to deal with. Therefore we 
use a different approach. We will use a reflection of r„ with respect to its main anti-diagonal, 
which turns the above ratios into expressions of the form 



det(z/„ 



det(2;/„ - Tn) 



Pn{z) 



Note that the matrix in the numerator of the left hand side, is obtained from the matrix in the 
denominator by skipping its last j rows and columns (rather than the first j rows and columns, 
thanks to the anti-diagonal reflection.) The limiting ratios for n oo will then be obtained with 
the help of equation (jl.2ip (depending on the residue class of n modulo p) . 

Now we work out the above ideas in detail. For a fixed /c e {0, . . . ,p — 1}, let 



T 



(fc) rp[l,...,k;l,...,k] 



i.e., T^^^ is the infinite matrix obtained by deleting the first k rows and columns of T . Of course, 
we understand that T'^'^^ = T. Note that T^*^' is also a tridiagonal block Toeplitz matrix 



V 









^(fc) 






Bp) 


Bq 


Bi'l 






^0 
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where the blocks of T^*^^ are also of size p x p and can be computed easily. Let Hp be the p x p 
permutation matrix that consists of I's in the main anti-diagonal and O's elsewhere, i.e. 



np(i,j) 



ii i + j — p 
otherwise, 



1, 



(3.1) 



where Ylp{i, j) represents the entry in row i and column j of Hp. We construct now a new block 
Toeplitz matrix T^*^^ as follows: 



f{k) 









Bq 








Bq 


Bi'l 






^0 



where 



B. 



(k) 



-1,0,1, 



where the superscript ^ stands for the matrix transposition. So the block B\ ' is obtained by 
refiecting B^ with respect to its main anti-diagonal. 
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(3.3) 



Lemma 3.1. Fix n ~ mp + k, with m G Z>o and k € {0, 1, . . . ,p — 1}. Then the Chebyshev- 
Nikishin polynomials {Pn-jYj^Q are obtained by the formulas 

Pn-j{z) = det(z/„_, - (f < J- < n, (3.2) 

where In-j is the identity matrix of size n — j and is the truncation of T^'^^ to the first n 
columns and rows. 

Proof. Formula (j3.2p follows directly from the relation T^'^ = n„ T„ 11 J , where n„ is the n x n 
permutation matrix defined as in p.ip . □ 

We keep n = mp + k fixed in what follows. If we set A = T^^), then Lemma [01 asserts that 
the polynomials {Pn-j}^^o 

of the form d';P , cf. So (see Remark lOl) we know 

that 

/ Pn-l Pn-l + Ci,iPn-2 Pn-1 + Cl^p-lP„-2 + C2,p-lPri-3 + ' • ' + Cp^i^p^iPn-p 

v~p~' p ' ' ■ ' ' P 

^ ^ n ^ n ^ n 

is a Hermite-Pade approximant to the system of functions 

(Jk) (k) Jk) (k) Jk) (fe) Jk) (k) 

[9l t9i +Ci,i52 i9l +21,232 + C2,253 ^■■■t9i +Ci,p-192 + r Cp-l,p-l5^ ' j, 

with respect to the multi-index (|1.2I) and the functions gf^ are the resolvent functions (12. 9p 
associated with the operator T'*^). 

4 Proof of Proposition 11.41 

In this section we prove Proposition 11.41 on the ratio asymptotics of the Chebyshev-Nikishin 
polynomials P„ . In the next section we will use this to prove Theorems 11.21 and 11.31 

Recall the notations in Section [L^ Since the block Toeplitz matrix T in (|1.8p - (|1.10p has a 
lower Hessenberg shape, each set Pfc in (|1.19p is non-empty and is a finite union of analytic arcs, 
see [TOl Prop. 1.1 and Lemma 2.1]. From P or [TOj Example 5.3] (or straightforward verification) 
we also know that the roots Zk{x) in (|1.18p have the asymptotics 



zi{x) ^ + 0[-p^], X ^ oo, 

^ ' (4.1) 

Ci|a;| < |zfe(a;)| < C2|x|, x ^ oo, 2<k<p+l, 

for certain constants Ci,C2 > 0. Therefore Fi is compact. 

We also need some results and notations from ISj . Let TZ denote the compact Riemann surface 



p 

'^k 

k=0 



Ti=\Jn 



formed by gluing in the usual crosswise manner the consecutive sheets 

7^o:=C\Al, 7^fc =C\(AfcUAfe+i), fc = l,...,p-l, TZp:=C\Ap. 

Let ^p^'^^ k = 1, . . . ,p, denote the meromorphic function on TZ whose divisor consists of a simple 
pole at oo^'^-' e TZq and a simple zero at oo^'^^ e TZk, with the normalization 

^P^''\z)=z + 0{1), z~^oo^°l 
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Here cxo*^*'^ denotes the point at infinity on the sheet TZk- Let ipj'^'' denote the restriction of V''-'^-' 
to the sheet TZj . The main result of [S] asserts that 

lim ^^^^^^^±^=:F,(.) = ^f (z), zeC\Ai, fc = l,...,p, (4.2) 

recall (jl.6p . From Liouville's Theorem we deduce that for any pair of indices j, k, the function 
is constant on TZ. In particular, Fj — Fk is also constant. 

From the recurrence 

zQunp+k Qnip+k-\-l ^" ^mp+kjnp+kQmp+k 4" ' ' ' ^ O^mp+k.mp+k—pQmp-'tk—pj 

dividing by Qmp+k and taking the limit as m oo, we deduce that for each fee {0, . . . ,p — 1}, 
z = J^fc+i + bk.k + h.k-i^ + bk.k^2-^r-k 1 'r hk.k^p— %, , (4.3) 

i'k ^k-fk-1 ^k- •■-fk-p+1 

with the obvious identifications hk,k-i = bp+k,p+k-i 

if z > fc and Fk-i = Fp+k~i iH > k (cf. (ITT^ 
and (II. 7p ). Multiplying both sides of (j4.3p by Ffe • • • F^-p+i, we obtain 

zFk ■ ■ ■ Fk-p+i — Fk+iFk ■ ■ ■ Fk-p+i + bk^kFk • ■ ■ Pk-p+i + ■ • ■ + h^M-p+iFk+i + bk,k-p- (4.4) 

Except for the term bk.k-p, all the other terms in (|4.4p contain a factor Fk+i = Fk-p+i. Writing 
all the functions j ^ fc + 1, as Fk+i plus some constant, we arrive at an algebraic equation 
of degree p + 1 satisfied by Fk+i{z) — xp^'^^^^z) of the form 

Fk+i {zf+^ + Cp{z) Fk+i {zf + . . . + Ci (z) Fk+i (z) + bkM-p - 0. 

By analytic continuation, this equation is also satisfied by the other branches V'^'^^^'' j * = Ij ■ • ■ jP- 
It follows that 

bk,k-p^{-ir+'fl4'^'\z), zee. (4.5) 

1=0 

In particular we have that fp ^ 0, see (jl.23l) . This allows us to write (ll.22p . 

Recall the strong asymptotics in p.24p . We remark that the function detF^P''''^^^{zi{x),x) 
in (|1.24p has at most a finite number of zeros outside Fi. Indeed, the contrary would mean that 
det fI^'-'^'+^I (zi(x), x) is identically zero as a function of x, which by analytic continuation would 
imply that each of the coefficients in Widom's formula (|1.22p for Pmp+k(x) is identically zero. 
But then Pmp+fc(a;) = 0, which is clearly a contradiction; see also [TUl Lemma 5.7]. 

From (|1.24l) we deduce that 

^^^^,=^d^) (4-6) 

and 

li^ ^^P+^(^) ^(_l). detFb^^+^l(z, (.),.) 

m™oo P^j,{x) ^ deiF\P^\z^{x),x) ' '•■•'^ ' ^ ' 

uniformly on each compact subset of C \ (Fi U .4) with A a finite set. 

We want to relate the ratio asymptotics in (14.71) and (|4.2p . Recall that all the zeros of all the 
polynomials Qn are contained in the interval Ai, and that the zeros of Qn and Qn+i interlace 
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(see [SJ Theorem 2.1]). We can then apply [7( Section 4] (taking into account that we do not 
know yet that Fi C M) and obtain for all a; G C sufficiently large that 



n— >oo 



Moreover, [71 Remark 4.2] shows that the vector 

Q(x) := lim 




satisfies the matrix- vector product relation 

F(zi(x),x)Q(a;) =0, 

for all a; S C sufficiently large. From the interlacing of zeros we know that all the limiting ratios 
VaiLm^oo Qmp+k{x) I Qmp(x) are analytic in C \ Ai. Thus we get by analytic continuation that 

lim-^^ = zi(x), a;GC\Ai, (4.8) 

and 

F{~zi{x),x)q,{x) a;GC\Ai, (4.9) 

where zi{x) is the (unique) root of the algebraic equation f{z,x) — that satisfies zi{x) ^ x~p 
for a; — > CX3 and that depends analytically on x G C \ Ai. Obviously zi{x) — zi{x) for all x large 
enough. We will see further that zi{x) = Zi{x) holds for all x £ C \ Ai. 

From the above discussion we see in particular that the root zi{x) can be analytically con- 
tinued from a neighborhood of infinity to C \ Ai. 

Applying Cramer's rule to (|4.9p gives us 

Qn.p+.(x) ^ , detFiP^>^^Xz,ix),x) 
J™oc Q„p(x) ^ ' detFb.il(zi(x),a;) ' l^.iUj 

for all fc = 0, . . . ,p — 1. Comparing this with (I4.6p - (|4.7|) . we see that the proof of (|1.2ip will be 
obtained once we know that zi(x) — zi{x) for all x G C \ Ai. 

Incidentally, we point out that an alternative proof of (j4.8p and (I4.10p could be constructed 
with the help of the generalized Poincare theorem in [23] , rather than using ^ . 

From and we see that 

'^^"^ = ^«(x).'.4^)(x)=^i(x).'.W ^'-''^ 

The right hand side can be analytically extended to a meromorphic function on the Riemann 
surface TZ, with restriction to the sheet TZk-i given by 

ik{x) := (1) , / .M , xGC\(Afe_iUAfe), (4.12) 

V'fe_i(a;)---V'fe_iW 

for fc = 1, . . . , P+ 1, where we understand that Aq = = Ap+i. Thus z\^(x) and Zfe+i(x) are each 
others analytic continuation across A^. From the uniqueness of analytic continuation we then 
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get that for each a; G C, 2i(x), . . . , Zp+i{x) are a permutatfon of the roots 2:1(0;), . . . , 2p+i(x) of 
the algebraic equation f{z,x) = 0. If we would be able to prove that 

\~zi{x)\<\h{x)\<--- <\2p+i{x)l xeC, (4.13) 

and 

\zk{x)\ < %+i{x)l a;eC\Afe, (4.14) 
then we could conclude from (|1.18p that Zk{x) — Zk{x) for all x G C \ (Afc_i U A^). 

The rest of the proof is devoted to proving (|4.13p - (j4.14p . We will do this by a rather intricate 
argument involving the total masses of certain equilibrium measures. An alternative proof may 
be found in |1] Propositions 2.1-2.2]. 

To start, we define a measure Sfe on A^ with density 

dsk{x) = — {^^^-^^^Adx, xeAk, (4.15) 
27ri \zk,+ (x) Zk-[x)J 

A: = 1, . . . where the prime denotes the derivative with respect to x, and where the + and — 
subscripts stand for the boundary values obtained from the upper or lower half of the complex 
plane respectively. Note that the density (I4.15P is well-defined in the interior of A^ , and blows up 
at worst like an inverse square root at the endpoints of A^ (this follows from the representation 
of Zk in local coordinates near the endpoints). We claim that Sk is a real- valued (possibly signed) 
measure on A^ with total mass 

Sfc(Afe) / dsk{x)=p+l-k, fc = l,...,p. (4.16) 
To prove Equation (|4.16l) . we first derive the following relation for the Cauchy transforms: 



L 



dsk-i{t)- -dsk{t) = ^^, xeC\{Ak-iUAk), (4.17) 



X -t X - t Zk[X) 



for k — 1, . . . ,p, where we understand that Aq = and sq — 0. Indeed, equation (I4.17P follows 
easily by contour deformation and by using that Zk-i.±{x) — Zk,^{x) for x G Afc_i. 
From 

zi{x) ^ x~P, X 00, 

Skix) ^ CkX, X — ?> 00, 2 < fc < p + 1, 

we get for the logarithmic derivatives that 

z'kjx) / -p/x, k = l, 



Zkix) \ 1/a;, k = 2,...,p, 

for X 00. Thus by equating the 1/x terms in the asymptotics for a; — > 00 in both sides of 
(|iT7)) we obtain 

-p, fc = 1, 



Sfc.i(A,_i)-«fc(A.) = | ^ fc = 2,...,p. 

The claim (|4.16p follows from this by means of an upward induction on fc = 1, . . . ,p. 
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The symmetry under complex conjugation shows that Zk{x) = Zk{x) where the bar denotes 
the complex conjugation. We then obtain from (|4.15l) - (|4.16p that 

1 /■ / z' (x)\ 

' Im J''^ ) dx^p+l~k, (4.18) 



with Im denoting the imaginary part of a complex number. 

On the other hand, consider the (positive) measure on Tk with density |10j 

1 /z'i.(x) z' (a;)\ 

This measure is well-defined on each open analytic arc of , recalling that Tk is a finite union 
of analytic arcs. Note that we are using the roots Zj and not Zj. The measure ak, restricted to 
the real line, takes the form 

1 / z' (x)\ 

dak(x) = - V Im ] dx, xeTkn R. (4.20) 

Now, let us take a fixed open interval J C R that does not contain any intersection points 
or endpoints of the analytic arcs constituting Ffc, for every k. We also ask J not to contain 
isolated intersection points of the sets with the real axis. Thus there exists an open connected 
set J7 C C such that J7 n R = J and moreover J7 n Ffc is either empty or equal to J, for any 
k = 1, . . . ,p. The boundary values Zk.+{x) for x e J are then uniquely defined and they vary 
analytically with x. 

On the interval J, there exist indices 1 < mi < m2 < . . . < tol < P such that 

|zi^+(a;)| = . . . = |z„jj,+(x)| < \z,n^+i^+(x)\ = ... = \z„r2^+{x)\ < ... 

< \zm^+i^+{x)\ = ... = |zp+i^+(a;)|, (4.21) 

for all X € J. We define toq ■= and mi+i := p + 1. 

We will see later that nik+i — ink G {1, 2} for all k, i.e., each "cluster" \zmk+i.+ {x)\ — . . . = 
\^mk+i,+ {x)\ in (I4.2ip can only have length 1 or 2. 

From the ordering (ll.lSp we find for each of the clusters \zmk+i.+ {x)\ = . . . = \zjn^_^-^^+{x)\ in 
(|4:2T|) that 



(Re logz„i^+i(x + ie)) 



< . . . < 

e=0+ 



(Re log Zm^^j {x + ie)) 



, X G J, 

e=0+ 



for any k — 0, . . . , L. Using the Cauchy-Riemann equations this implies 



Inr ff™^±i44^ >--->Im f^^^^^^^V ^ ^ J, (4.22) 

k = 0,...,L. 

By the symmetry with respect to complex conjugation we have a pairing of the numbers in 
(|4.22p : they appear in positive-negative pairs in the sense that 

Im '—^ = -Im i — , J = 1, . . . , mfe+i - m^, (4.23) 

\^rnk+],+ \X) J \Zjn^ + i + l-j,+ (X) I 



14 



for any a; G J and k = 0, . . . , L. Note that the numbers in (|4.22p can also be identicahy zero, in 
the case of a real- valued root. Clearly there is at least one real root if mk+i — ruk is odd. 



Let TOfc : 
cTThk in that 



2 



From the above considerations we find for the density of the measure 



dx 



— > Im 

TT ^ 

- F 



i=mfc + l 



1,+ 

I 

z 



[x) 
[x) 



X e J, 



for any k — 0, . . . , where the first equality follows from the cancelations arising from (j4.23p . 
and the second equality follows from (j4.22p - (j4.23l) . Hence the sum of the total masses of all the 
measures Uk over the interval J can be bounded from below by 



cri{J) + . . . + (jp{J) > a,n„{J) + ...+ 

P+i 



1 
2^ 



Im 



Zj,+ {x) 



(4.24) 



dx. 



Putting A := IJ^^^ we then obtain 

ai(A) + ... + ap(A) > ^ I 

7a 



Im 



-3.+ 



[x) 



[x) 



dx 



(4.25) 



On the other hand, the (positive) measure Cfc satisfies [10, Cor. 4.2] (or see the proof of (|4.16l) ) 



crfc(A) < crA;(C) = crfc(rfe) =p 

for k — 1, . . . ,p. Using this in (|4.25l) we get 



1 - fc, 



(4.26) 



Y^ip+l-k) > 



k=l 



> 



E<^fe(A) 

fc=l 

1 

2^ 



A 



1 

2^ 



A 



p+1 

E 

k=l 

p+1 

E 

fe=i 



> 



-t 

k=l 

p 



Im 



Im 



Im 



fc). 



-k,+ 



(x) 



Zk,+ ix) 
Zk,+ ix) 



(x) 



Zk,+{x) 



dx 



dx 



dx 



k=l 



where the second relation is (I4.25P , and the third one follows since the roots Zk form a permutation 
of the roots Zk- The fourth relation uses that on A^ there are precisely two non-real roots Zk{x) 
and Zk+i{x), which are complex conjugated. Finally the last inequality is (j4.18p . 
From the above chain of inequalities we find: 
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Lemma 4.1. (a) We have 



p 

rfeCA=|jA,, fc = l,...,p. (4.27) 

(&) Clusters of length > 3 m (|4.2ip cannot occur unless they have all entries in (|4.22p equal to 
zero. 

Proof of Lemma \4-.1\ (a) By comparing the two outermost terms in the above chain of inequal- 
ities, we see that each of the inequahties > is actuaUy an equahty. Since the first inequaUty in 
the chain comes from (I4.26p . this imphes in particular that crfc(A) = CTfc(C) = akij'k) = p-\-l — k 
for all k. This implies (I4.27p . since each analytic arc of Tk has a strictly positive mass under at- 
Let us prove the statement made in the last sentence. We note that the product zi{x) . . . Zk{x) 
is analytic in C\rfc |13[ Prop. 3.5] and moreover, by splitting the analytic arcs of Tk into finitely 
many sub-arcs if necessary, the analytic continuation of zi (x) . ■ . Zk (x) across each such sub-arc 
is of the form Zj^{x) . . . Zj^{x) with {ji, . . . ,jk} ^ {1, . . . ,k} [l3l Prop. 3.4]. This then implies 
that log \ zi{x) . . . Zk{x)\ is nowhere harmonic on F^. Since — log \zi{x) . . . Zk{x) \ + C (for some 
constant C) is the logarithmic potential of the measure crfc in (j4.19p ^TU[ Prop. 5.10], we then 
conclude from the non-harmonicity that ak is indeed nontrivial on each sub-arc of F^. This 
proves the claim. 

(b) From the strictness of the inequalities > in the above chain of inequalities, we also get 
that the inequality > in (j4.24l) is an equality. This is easily seen to imply Part (b). □ 

Now we are in a position to prove (|1.20p . First we will prove that Fi = Ai. From (14.271) we 
already know that Fi C A. Then the equality zi{x) = ^1(2;), which is known to hold for all x 
sufficiently large, remains valid for all a; G C \ A. Using this, we find 

zi^+{x) = zi.,+ix) = zi,-(x) = 2:1, _(x), a; e A \ Ai, 

where we used that zi(x) is analytic away from Ai. Hence zi{x) is analytic on A \ Ai and then 
[HI Prop. 3.4] implies that (A \ Ai) n Fi = a nd so Fi C Ai. The reverse inclusion Ai C Fi is 
immediate since zi^^(x) = z\^j^(x) — z\^-{x) = z\^^{x) for x G Ai, yielding two different roots 
with the same modulus. This proves that Fi = Ai. 

From the non-triviality of a\ on Fi = Ai (see the proof of Lemma [4.1f a')'). we have that 

Im f fi^) J, for almost aU a; G Ai. Thus the cluster Im { fLt^f)) > . . . > Im f ^^^^i^) on 

the interval Ai has length mi < 2, in view of Lemma |43Ib). On the other hand, we also have 
that TOi > 2, since z\^±(x) = Z2,ip(x) are complex conjugate for x G Ai. Thus m\ = 2. This also 
shows that Z2{x) = Z2{x) for all x near Ai and therefore for all x G C \ A, since F2 C A (again 
from (|T?7)) V 

Now we prove that F2 — A2. In view of the last paragraph we have 

zi,+ (x)z2,+ (a;) = zi^+{x)z2,+ {x) = zi _(x)z2,_(x) = zi _(a;)z2 _(a;), a; G A \ A2, 

where we used that the product Z1Z2 is analytic away from A2. Thus Z1Z2 is analytic on A \ A2 
and so 13, Prop. 3.4] implies again that (A \ A2) fl F2 = and therefore F2 C A2. The reverse 
inclusion A2 C F2 follows again since Z2.+ (x) — Z2^+{x) — Z2,-{x) = Z2^-{x) for x G A2, yielding 
two different roots with the same modulus. This proves that F2 = A2. 

This reasoning can be extended to show that F^ — A^ for each fc = 3, . . . ,p. In particular, 
we have that Zk{x) = Zk{x) for x G C \ (A/c_i U A^). This concludes the proof of Prop. 11.41 □ 
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5 Proof of Theorem 11.21 

Using (|1.21l) and a telescoping product we deduce that 
1- Pmp+k-j (x) _ 1 

where each Fj is understood as -Fjmodp- It follows that for n = mp + k, < k < p~l, the vector 
p.3p is a Hermite-Pade approximant to the system of functions 

gw = (g(^-),g(''^...,gW), 

where 

In particular, we obtain the following orthogonality conditions 

Pnix) x'Gf \x)dx^O, l = 0,...,n,-l, j = l,...,p, (5.2) 

where 7 is an arbitrary closed Jordan curve surrounding Ai. From (j5.2l) we see that the polyno- 
mials Pn satisfy certain multiple orthogonality relations, but with the weight functions G^'^'' (x) 
depending on the residue class k oi n modulo p. To get around this issue, we will make a clever 
choice of the constants Cij^i in (|5.ip . 

Proposition 5.1. Let j and k be fixed. The constants c.i.j_i in (IS.ip can be chosen so that 

Gf\x)^ \ (5.3) 
^ Fk-j+i{x) 

where again Fk-j+i is understood as F(^i^^j^i^ modp- 

Proof. The function Fk{x) in (jl.6l) has a simple pole at x = 00, with 

Kk := lim Fk{x) - x = -bk-i.k-i, (5.4) 

a;— >oo 

by virtue of (|4.3p (with z replaced by x). As already mentioned, Liouville's theorem implies that 
ipU) — ^(k) constant on TZ, for any j, fc G {1, . ■ . ,p}. Combined with (|5.4p and (14. 2p this yields 

F,(x) - Fkix) = V4'^(a;) - ^^''^(x) =Kj~Kk, xeC\Ai. (5.5) 

Now let j and fc be fixed. We will evaluate the sum (|5.ip from right to left. Formally, we do 
this by defining functions yi{x), I = 0, . . . , j with yo{x) — and then recursively by 

yi{x)=yi-i{x)+cj-i,j-i—— --, (5.6) 

J'kyxj-t'k-iKX) ■ ■ ■ Pk-j+i[x) 

for I = 1, . . . , j , with cqj-i := 1. So yi{x) is the sum of the last / terms in (15. ip . and in particular 
y,ix)^Gf\x). 

We will prove by induction on I that the constants Cj_;.j_i can be chosen so that 

?/;(a;) = Cj_;.j_i— — — -— --, (5.7) 

Fk{x)Fk-i{x) ■ ■ ■ Fk-j+i+i{x)Fk-j+i{x) 



17 



for any I = 1, . . . ^ j . By taking I = j this proves the proposition. 

It remains to prove (|5.7p . For I = 1 this is trivial. Now assume by induction that (|5.7p holds 
with I replaced by I — 1 G {1, . . . ,j — 1}. From the definition (15.61) we then get 

1 1 

yi(x) - c,-,+i^,_i . . . F^^^^^^^^Fk-,+iix) ^ '"'-'^'-^ Fk{x)Fk^i{x) ■ ■ ■ Fk^,+i{x) ' 

Putting on a common denominator yields 



Cj-i+i,j-i + Cj-ij-iFfc_j+i(x) 



Fk{x)Fk^i{x) ■ ■ ■ Fk-j+i{x)Fk^j+i{x) ' 
Now by taking Cj_;+i.j_i Cj^ij^i{Kk-j+i — Kk-j+i) and using (jS.Sp . this becomes 



yi{x) = Cj-;j_ 



'-'^'-^ Fk[x)Fk-i{x) ■ ■ ■ Fk-j+i{x)Fk-j+i{x) ' 

which shows that (j5.7p holds for /. This proves the induction step. □ 

From (j5.2p and (|5.3p . we see that the Chebyshev-Nikishin polynomials P„(x) satisfy the 
orthogonality conditions 

Pn{x)x' dx = 0, l = 0,...,nj -1, j^l,...,p, (5.8) 
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with again 7 an arbitrary closed Jordan curve surrounding Ai. Indeed, if n = mp + k with 
£ {0,l,...,p — 1}, then (15.81) can be checked for each fixed residue class k individually, using 
the above considerations and the fact that the orthogonality measures with the same number of 
orthogonality constraints can be freely interchanged. 

By shrinking the contour 7 in (j5.8l) to the interval Ai and using the Stieltjes-Perron inversion 
formula, we get (see (jl.l4p ) 

P„{x) di^j{x) =0, l = 0,...,nj-l, j = l,...,p. (5.9) 

Ai 



This proves Theorem 11.21 □ 



6 Proof of Theorem 11.3 



The proof uses some ideas from Aptekarev-Kaliaguine-Saff [3], see also [TT| . 
Write (z/f\...,4^^) := 

(i^i, . . . , I'p). These measures will form the first layer of the Nikishin 
hierarchy. Denote the Cauchy transform of i^j^"* by 

F^'\x):=l/F,{x), J = l,...,p, 

recall (|1.14p . By the Stieltjes-Perron inversion principle, the density of Vj^^ is given by 

'-(4^(.)~4^(x)), xeA,, 



dx 2Tii 
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with again ± denoting the boundary values from the upper and lower half-plane respectively. 
With the help of (14.21) this becomes 

where the last equality uses that, according to (j5.5L 

^^'\x)::=^tl)''^\x)+Kj-Kk, l<j,k<p, (6.1) 

which implies 



(Notice that ^/'o'^Vl''^ € 'H(C \ As), i.e., it is holomorphic in C \ A2.) 
Let us prove that the ratio of densities of z^j^'' and 

Fl^{x) - Fl[{x) 

can be analytically extended to C \ A2 and can be written as a Cauchy transform 

f ^di.f(i), :^eC\A2, j = 2,...,p, 

(2) (2) {2^ 

for some measure with constant sign supported on A2. The measures (^ ? • ■ • 7 '^p ) form 

the second layer of the Nikishin hierarchy. 

In fact 

^ ^ ^"4^)(x)V^«(x)' ^ 

and thus 

F,^^' = |y^eH(C\A2), ,=2,..., p. (6.2) 

Moreover, at infinity the numerator takes a finite value whereas the denominator has a simple 
pole, therefore 

Ff^(x)=oQ^, x^oo, 

having a simple zero at 00 and no other zeros in C \ A2 . 

Let 7 be a positively oriented closed Jordan curve surrounding A2. By Cauchy 's integral 
formula we have that 



for all X exterior to 7. Shrinking 7 to A2 we get that 



F^'\x) 

where 



X — t 
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Since ^/'^''^ G ^(0 \ Ai), fc = 1, . . . ,p, and Ai n A2 = 0, the density of vf' takes the form 



dx 2m [^^'-^''^ ^^'^^''n 27rz^O-)(^) i^},(x)) ' ^ 

On the other hand, using (j6.ip . the expression in the last parentheses reduces to 

since = 4'2+ix)'4'i\.{x) = V-'i"'-(a^)V'2"'-(^) can be extended analytically to a 

neighborhood of A2. Consequently, 

da; 2m 4^\x)4^\x)4^\x) ' ^ 

where V'o"' Vi"' ^2"''' 'H(C \ A3) and £ H(C \ Ai). Moreover, from the symmetry of these 
functions we deduce that ip'^^ / ip^^ takes real values with a constant sign on A2. On the 

other hand, 



4^1 (x) - 4iXi^) = i'llix) - 4i -i^) = 2ilm{tl;[^l{x)}, X e A2. 

Should 4il{x) = for some X in the interior of A2 , we would have that at that point (^) ~ 

4i +{x). This is clearly impossible because V'*'^' is one to one on TZ whereas x+ and x^ are 

(2) 

distinct points on this Riemann surface. Consequently for each j — 2, . . . ,p, the measure Vj 
has a constant sign on A2. 

In general, fix I e {2, . . . ,p} and assume that we have defined measures (t'^^j^^-', . . . ,i>p ^■') 
supported on A;_i, forming the — 1) layer of the Nikishin hierarchy. Denote the Cauchy 
transform of by j ~ I — 1, . . . ,p. Assume that we have also shown that 



mi-2,i-2 ip''^>{x)...ipY-2i^) 



and 



mi-ij 4o~^^ jx)... ipfZ^^ {x) ^^,^(1) ^^^(1) 



{i>\'\Ax)-4'\4x)), 



for certain constants mi_2j, "t-z-Ij" £ K \ {0}. (In fact, toqj = Ij J = 0, mi.j 
j = 1, . . . ,p, and TO2.J = - Ki,j ^2,... ,p.) 
Then, we set 

(x) ■= ^TJ—T) TT-n ' xeAi_i, 

Ft,lix)~Fl[-'lix) 
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for j = I, . . . ,p. We will see that this function can be analytically extended to C \ Aj, and that 
it is again a Cauchy transform 

^f(^)=/ ^di/f(t), xeC\A,, (6.3) 

for some constant sign measure v'j'^ supported on A;. The measures ■ ■ ■ , Vp'^) form the Zth 
layer of the Nikishin hierarchy. With this we will conclude the induction. 
In fact, using the formulas from the induction hypothesis, it follows that 

Fy(x)^ -T -T , a; e A/_i, 

' m-u-i 4^\x)...4'l'M^) 

which by the properties of the branches 'ipk~^^ and ^\^\k — 0, — 1, can be extended 
analytically to all C \ A; . At infinity the numerator takes a finite value whereas the denominator 
has a simple pole. Consequently, 



oo. 



(x) e7^(C\A0, F]'>{x)=0\^- 
Therefore, using the Cauchy integral theorem, we obtain 

Im X — t 

for all X exterior to 7, where 7 is a positively oriented closed Jordan curve that surrounds A;. 
Shrinking 7 to A; , we find that (|6.3D takes place with 

d'^f (^) 1 fM(\ M(\\ V'o"'^(^)---#V^(^) ^(1) .X ,(1) .A 



(6.4) 

for X € Ai, with mi_j — -^'-1) ^ Here, we also used (j6.ip . Using the expression for 



— ^j^;^ — , arguing as we did above for the case I — 2 one deduces that ^'j has constant sign on A; 
for each j — I, . . . ,p. The fact that the measures are absolutely continuous with respect to the 
Lebesgue measure is a direct consequence of the expression we have found for them. This ends 
the proof of Theorem 11.31 □ 



Remark 6.1. Note that the generating measures pj in Theorem ll.3l are nothing but the "diagonal" 
measures in the Nikishin hierarchy, i.e., pj — given by (16.41) . 

Remark 6.2. The recursive relation given above between the constants mij allows us to write 
for ain = 2, . . . ,p and j — I, . . . ,p, where the second equality follows from (|5.4p . 
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7 Proof of Theorem 11.11 

Let 

OO 

Fk{z) = z - bk-iM-i +'^—, |z| > ^0, k = l,...,p, 
1=1 ^ 

denote the Laurent expansion at infinity of Fk, where we are using (j5.4p - (|5.5p . Set 

OO 

1=1 

Obviously, 

Fkiz)^ Fiz)-bk-i,k-i- (7.1) 

Since -F is a translation of Fk , which is one branch of a conformal representation of the Riemann 
surface TZ onto C (see (14. 2I) '). F itself extends to a conformal representation of TZ onto C. We 
denote by the analytic extension of F to all TZ. For j ^ k we have that ^ bk-i,k-i, 

otherwise Fj — Fk. Taking these facts into account, we conclude that there exist p distinct points 
Co, ... , Cp~i £ such that 

-^(Cfc) = 6fc,fc, fc = 0, 

In the sequel C„ = if n = A; mod p, fc G {0, . . . ,p — 1}. 
It follows from and ((7?T|) that 

bk,k-i = lim {z - F{z))Fk{z) = -ci =: /3, /c = 0, . . . ,p - 1, 

so (jl.lip is proved. 

Using again (|4.3p and (|7.ip . we find that 

zFfc • ■ • Fk-pj^i — ■ ■ ■ Fk-p+i + bk,k-iFk-i ■ ■ ■ Fk-p+i + • ■ • + bk,k-p- (7.2) 

Writing equation (|7.2p substituting fc by fc + 1 and taking into consideration that Fk ■ ■ ■ Fk-p+i = 
Fk+i ■ • ■ Fk-p+2, we obtain the equality 

bk,k-iFk-i ■ ■ ■ Fk-p+i + • • • + bk^k-p = bk+i,kFk ■ ■ ■ Fk^p+2 + • ■ • + bk+i^k-p+i, 
which due to (|7.ip , and taking into consideration that F extends to all TZ, takes the form 

bk,k-i{J' — fefe-2,fe-2) ■ ■ ■ {J^ — hk~pM-p) + • • • + bk^k-p+i{J' — bk-pM-p) + bk,k-p = (7.3) 

bk+i,k{^ — bk~i,k~i) ■ ■ ■ {J- ~ bk-p+iM-p+i) + ■ • • + bk+i,k-p+2{J' — bk-p+i,k-p+i) + bk+i,k-p+i- 
Evaluating (|7.3p at Cfc-p+i G TZ, it follows that 

bk,k-p+i{bk-p+i,k-p+i — bk~pM~p) + bk,k-p = &fc+i,fc-p+i, (7.4) 

which is equivalent to (jl.l2p when I — p ~ 1. 

Now, evaluating (|7.3I) at Ck-p+2 G TZ, we obtain 

bk,k-p+2{bk-p+2,k~p+2 — bk~p+l,k-p+l){bk-p+2,k-p+2 ^ bk-p,k-p) + 
bk,k-p+l{bk~p+2,k-p+2 — bk-p,k-p) + bk,k-p — 
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bk,k~p+2{bk-p+2,k-p+2 — bk~p+lM-p+l){bk-p+2,k-p+2 — bk-p,k-p) + 
bk,k-p+l{bk-p+2,k-p+2 — bk-p+lM-p+l) + 

bk,k-p+i{bk-p+i,k-p+i — bk-p,k-p) + bk,k-p = 

bk+l,k-p+2{bk-p+2,k-p+2 — bk-p+l,k-p+l) + bk+l,k-p+l- 

Using (|7.4I) . this equality reduces to 

bk,k-p+2{bk-p+2,k-p+2 — bk-p+l,k-p+l){bk-p+2,k-p+2 — bk-p,k-p) + 
bk,k-p+l{bk-p+2,k-p+2 — bk-p+l,k-p+l) = bk+l,k-p+2{bk-p+2,k-p+2 — bk-p+l,k-p+l) ■ 

We can cancel out bk-p+2.k-p+2 — bk-p+i,k-p+i{¥' 0) from both sides and we obtain 

bk,k-p+2{bk-p+2,k-p+2 — bk-p,k-p) + bkM-p+l = bk+lM-p+2 C^-S) 

which is prr^ for ? = p - 2. 

In general, assume that (I1.12p is valid for all I ^ {p — l,p — 2, ... ,p — m}, with m < p — 2. 
In order to prove that it is also valid for I p — m — 1, we evaluate (j7.3l) at C/j^p cind we get 

^k,k~-T(^k~-Lk-T^ ^k-T-l,k-I-l) ■ ' ' (^k-T,k-T^ bk-p,k-p) + 
^k,k-li^k-l,k-T~ ^k-T-2,k-T-2^ ' ' ' ^^k-T.k-T^ bk-p.k-p) + 

■ ■ ■ + bk,k-p+2{bf._jf._j— bk-p+i.k-p+i){bj._jj._j~ bk-p,k-p)+ 

bk,k-p+i{bf^_jj._j ~ bk-pM-p) + bk,k~p — 
^k+i,k-ji^k-T,k-T^ ^k-T-i,k-T-i) ■ ■ ■ i^k-T.k-J- bk-p+i,k-p+i)+ 
^k+i,k^T~ii^k~T,k-T^ ^k-T-2,k-l-2) ■ ■ ■ i^k-lk-T'' bk-p+i,k-p+i)+ 

■ ■ ■ + bk+l,k-p+3{b}^_Ji^_J— bk-p+2,k-p+2){bf._Jj,_J— bk-p+lM-p+l) + 
bk+l,k-p+2{bf._J^._J— bk-p+l,k-p+l) + bk+l,k-p+l- 

On the left-hand side we replace bk,k-p+iibf._jj._j - bk-p,k~p) by 

bk,k-p+i{bf._j,^_i— bk~p+i^k-p+i) + bk,k-p+i{bk-p+i,k-p+i ~ bk-p^k-p)- (7.6) 

Now we apply (|1.12p for / = p — 1 , and this allows us to delete the second term in (17. 6p along with 
bk,k-p on the left-hand side and bk+i.k-p+i on the right-hand side. Now the factor bj^_j j ~ 
bk-p+i^k-p+i appears in all the terms on both sides of the resulting equation. We cancel out this 
factor everywhere to obtain a new equation in which on the left-hand side we have the expression 

bk,k-p+2{b]^_Y^]^_l — bk-p^k-p) + ^fe,fe-p+i, 
and on the right-hand side we have 

bk+l,k-p+z{bf.Jij._J~ bk-p+2,k-p+2) + bk+l,k-p+2- 

In a similar way, we now replace bk,k-p+2{bf._jf^_j- bk-p,k-p) by 

bk,k-p+2{bf._Jf._J— bk-p+2,k-p+2) + bk^k-p+2(bk-p+2,k-p+2 ~ bk-p^k-p) (''•7) 

and apply (|1.12p for I — 2. This allows us to delete the second term in (|7.7p along with bk,k-p+i 
on the left-hand side and bk+i,k-p+2 on the right-hand side. The factor b^_j^,_j— bk-p+2,k-p+2 
appears in all the terms on both sides of the resulting equation. We cancel out this factor. 

It is clear that continuing in this fashion we will arrive at the desired equation. Relation 
(|l.lip may be regarded a special case of (|1.12p for Z = 0. We are done. □ 
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Question 7.1. Suppose that a collection of real numbers {bk.j} is given satisfying (jl.lip - (ll.l2l) . 
What other conditions must be added in order that the given collection of numbers corresponds to 
the periodic coefficients of the recurrence relation satisfied by a sequence of Chebyshev-Nikishin 
polynomials? 

8 Formulas and strong asymptotics for the second kind 
functions 

In this section we obtain strong asymptotics and an exact Widom-type formula for the second 
kind functions 5'„^; in (|1.25p - (|1.26p . The function \['„_; (with 1 < I < p) satisfies the following 
boundary value problem: 

Boundary value problem 8.1. 

(a) '^n,i 'is analytic for z G C \ A; . 

(b) satisfies the jump relation ("^n.oiz) = Pn{z)) 

{^nA^})- - = 27rz vl/„ x £ A;. (8.1) 

with ni in (jl.2p . 

(d) ^„.i stays bounded near the endpoints ai and /3i of the interval A;, i.e., '^n.i{z) = 0{1) as 
z ^ ai and z ^ Pi. 

Properties (a) and (6) follow immediately from the definition of Property (c) is a 

consequence of the orthogonality properties satisfied by ^n.i-i with respect to the measure dp;, 
proved in [TS] Proposition 1]. Finally, [d] follows from the fact that "^n.i-i is analytic on A; and 

dp;(x) _ 'iyf \x) { 0{\x- a;|i/2), X ^ a;, 
dx " dx "I 0{\x - x^Pi, 

see (|6.4p and Remark 16.11 Given ^„,;_i (see (|8.ip ). it is readily seen that 'I'„^;(z) is uniquely 
determined by the above boundary value problem. 

Proposition 8.2. Fix I G {0, . . . ,p}. For any m large enough and for any fc = 0, . . . ,p — 1, the 

following Widom-type formula holds: 

(-l)P+fc ^ dctF[P-*^+il(zj(a;),a;) 
fp j=i+i ]Yttl,t^j{^ji^) - zt{x)) 

The constants ?Tii+i.i+i are given in (|6.4p - (|6.5D . 

Remark 8.3. Formula (j8.2l) can be computed for all a; G C \ ^, where A consists of those points 
a; G C for which Zi{x) = Zj{x) for certain i ^ j. Thus ^ is a finite set and it is formed by the 
zeros of the discriminant of /(z, a;) — in (11.171) with respect to the variable z. Note that the 
endpoints of the intervals A; are points in A. 



n--+^^-+i (0/,,(5irM (8.2) 
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The product of the entries rrii+i i+i in (|8.2I) can be evaluated using (|6.5p : 

l-l l-l 1-2 

1=0 i=l 1=0 

(The minus sign conies from mi.i = —1.) It also follows from (|1.2ip . (|1.24p . and (14.21) that 



^^-i^"^^" detFM(z,(a:),a:) ' ^^"^^ 



for all fc = 1 , . . . , p — 1 and j = and 



for all j = 1, . . . ,p + 1, so in fact all the quantities in (j8.2l) depend on the minors of the block 
Toeplitz symbol F(z,x) (|1.16p evaluated at the roots z = Zj{x). 

We will prove formula (j8.2p only for m large enough, although it may actually be valid for 
all m. 

Proof of Provosition \8.IA We prove (|8.2p using induction on I. The case I = reduces to (|1.22p . 
Let us assume as induction hypothesis that (|8.2I) holds with I replaced by I — 1. We need to 
check that the right hand side of (|8.2p satisfies the four conditions {a)-{d) in the boundary 
value problem 18.11 Actually we will prove that the right hand side of (|8.2p , which we denote by 
^mp+fc,i, satisfies a weaker version of (a)-(d) where these conditions are replaced by 

(a') 'i'mp+k,i{z) is analytic for z £ C \ (A; U A). 

{b') ^mp+fe,i satisfies the jump relation 

{^rnp+k,lix))- - ($mp+fc,i(x))+ = 27ri ^ mp+k,l-lix) ^^^^^ , X G A; \ A 

(c') As z — > oo, 

where (5 G R is some fixed constant independent of m. 

(d') Near each point G ^ there is a fixed integer qi such that ^mp+fc,i(-z) = — ai)"''/^). 

The integer will be precisely the multiplicity of a; = Oi G ^ as a zero of the discriminant of 
f{z,x) = 0. The solution to the boundary value problem {a')-{d') is unique only for m large 
enough, since then the (possible) poles at the points in A cannot compete with the high-order 
zero at infinity. This will prove that ^'mp+fc,; = ^mp+k,i, that is, (18. 2p holds. 

In view of (j8.3p - (|8.4p . the function ^mp+fc,i is a symmetric function of -zj+i, . . . , ^p+i, and it 
is also a symmetric function of the variables zi, . . . , zj. More specifically, if we view "^mp+k.i as 
a function A(zi, . . . , z;, z;+i, . . . , Zp+i, x) in the variables z^ and x, then 

A(zi, . . . , Z;, Zi+1, . . . , Zp+1, x) = A(zJ', . . . , Z;*, Zf^2: • • ■ 1 ^p+li 

where (zj, . . . , z^*) is any permutation of (zi, . . . , zj), and (z;*(_]^, ■ . ■ , -Zp+i) is any permutation of 
(zj+i, . . . , Zp+i). This property readily implies the analyticity of ^'mp+fc,; on C \ (A; U A), which 
is Part (a'). 
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For Part (b') we distinguish in the expression of '^mp+k,iix) the term with j = I + 1 from the 
terms with j ~ I + 2, . . . ,p + 1. We start with the latter terms. First we write the factor with 
i = / — 1 in the product in (j8.2p as 

W-i(^) - V'f-i(^) ^ W-i(^) " ^j-i(^) _ 1 l_ ,35. 

Observe now that the + and — boundary values on A; \ ^ of the expression 

are equal. Indeed, this follows since, with the identifications (j8.3l) - (j8.4l) . equation (j8.6l) depends 
symmetrically on the variables z/4.2, . • . , ^p+i, depends symmetrically also on zi, . . . , z;_i, and 
does not depend on zi or z;+i. 

The equality of the + and — boundary values of (|8.6I) also holds if we drop the factor 
l/'0j'2i(a;)- In view of (I8.5p . it then follows that if we subtract the + boundary value from the 
— boundary value of the terms in ^ mp+k.i{x) with j = / + 2, . . . ,p + 1, we obtain 

'^1,1 m m 

4'l^{x)^f\x) 

Here we also used that the product i^fli _{,x)iljf\ j^{x) — ■il^f}^^{x)-il^f\x) is analytic across A;. 

Using the relations Zi+i^±(a;) = z;.zp(x), il)\%{x) = '4>\-\^{x)^ ^ G and other simple 
considerations, it is easy to see that the difference of the — and + boundary values of the term 
in 'i/rnp+k.i{x) corresponding to j = Z + 1 is given by the expression in (j8.7p with the sum J2^td+2 
replaced by X^'^r conclusion. 



^j=l- 



4'\{x)4'\x) 



X Hj!i!g det^..^^^^i(z,(.),.) fn-.-..^C!^^^#^l ^.i^r-'- m 

By the induction hypothesis, we know that 'i'mp+k,i-i{x) is given by the expression in the last 
line of (|8.8p with tpfH^) replaced by V'f ~^^(a;). This observation and (|6.4p gives Part (6'). 

Finally, the conditions (c') (cf. (j4.ip ) and (d') are obvious. This proves (|8.2p . □ 

Remark 8.4. The functions defined by 
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where the measures vi are the orthogonahty measures in Theorem II. 2[ are also often cahed 
second kind ftmctions (see (THl pg- 672]), and represent the hnear remainders in the Hermite- 
Pade approximation to the functions J^^ '^z-t^ ' Z = 1, . . . ,p. It is easy to see that for aU m large 
enough, we have the formula 

^ ,,,M-tir^Y detF[^.^^^l(.,(x),a:) ^i^y-^ (8 9) 

We leave the justification of (|8.9I) to the reader. 

Corollary 8.5. The following strong asymptotic formulas hold uniformly on compact subsets of 
the indicated regions. For each fixed k G {0, . . . ,p — 1} and I G {1, . . . 

lim *™p+fc,;(a;)2:;+i(a;)'"+^ 

(-1)^+'- detfb.'-+^l(z,+i(x),x) ffl V.W(x)-^W(x) \ . 

nt=i,M/+i(^'+i(^)-^*(^)) \r=o 

and 



lim $„p+fej(x)z2(a;)™+^ 

_ (-1)^+^ detF[p^'=+il(z2(a;),x) 



/p neM/2(^2(^)-^t(x)) \^p^\x) ^«(a.). 



a;eC\(AiUA2). 



Proof. This is an immediate consequence of (|8.2p . (|8.9p . and the fact that |z;+i(a;)| < |zj(x)| for 
all j > ; + 2 and a; e C \ (A, U A/+i) (cf. (fLT5)) and ^^). □ 



9 Alternative systems of multi-indices 

The results in this paper were formulated for the standard "staircase" system of multi-indices 
n in (|1.2I) . More generally, let 11 := (tti, . . . , tt^) be an arbitrary permutation of (1, . . . ,p). For 
each n G Z>o, define 

(«?,...,„") ^^P^ 

as the unique multi-index such that 

<>... ><><-!, and In"! + + = (9.1) 

For the trivial permutation 11 = (tti, . . . , tt^) = (!,... ,p), this definition reduces to (|1.2I) . 

Let (5^(2) be the monic multiple orthogonal polynomial of degree n satisfying (|1.3I) . with 
Uk replaced by n^. The polynomials Q^{z) satisfy a recurrence relation of the form (jl.4p . 
with recurrence coefficients ajj^^^ having periodic limits as in (|1.7p . We define the block 
Toeplitz matrix and the Chebyshev-Nikishin polynomials P"(z) as in prSjl - pTTU)) and pTTO)) 
respectively, with instead of 6.;^ . 

We have the following generalization to Theorem 11.11 The numbers bfj satisfy the relations 
(|l.lip - (|1.12p . with (3 independent of the permutation IT. 

We also have the following generalization to Theorem 11.21 
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Theorem 9.1. The Chebyshev-Nikishin polynomials P„ satisfy the orthogonality conditions 




? = 0, . . . ,rifc - 1, k = l,...,p, 



(9.2) 



with Vk the measure defined by (jl.l4l) (independently of the permutation li), and with n]^ the 
k-th component of in (|9.ip . 

Theorem 19.11 can be proved in a similar way as Theorem 11.21 and we leave the details to the 
interested reader. See also the next paragraphs for more information. 

We define F^{z) as in (|1.6|) . with each Qn replaced by Q^. The main result in [5 implies 



th&t F^{z) = F^,{z). 

Define the quantities F^{z,x), zj^{x) and V]} as in prTC)) - prTg)) . with bfj instead 



of bi,j. Reproducing the proof in Section HI we are led to the following analogue of (14. lip : 



for all X sufficiently large, where we used the above observation that F^{z) = F.^^ (z). So for all x 
sufficiently large, z\^{x) = 2:1(2;) is independent of the permutation 11. Analytic continuation then 
implies that each of the roots zj}{x) = Zk{x) is independent of 11. This implies in turn that the 
algebraic equation f^{z, x) = f{z, x) = is independent of 11 and therefore also = = A^. 

Question 9.2. Given the intervals Ai, . . . , Ap, we have now obtained p\ block Toeplitz matrices 
T^, labeled by the permutations 11 of (l,...,p), which all have the same algebraic equation 
f{z,x) — 0, and which all satisfy (|1.20l) . Are there any other block Toeplitz matrices T of the 
form (jl.8l) - (jl.l0p (with arbitrary entries bij) leading to this same algebraic equation f{z, x) = 0? 
(or equivalently, for which Tk — A^ for all fc, using the notations in (ITTT6)) - ((ri9)) ?) 

We conjecture that the answer to the above question is negative. 
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